Phase of Ising spins on modular networks analogous to social polarization 
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Coordination processes in complex systems can be related to the problem of collective ordering 
in networks, many of which have modular organization. Investigating the order-disorder transition 
for Ising spins on modular random networks, corresponding to consensus formation in society, we 
observe two distinct phases: (i) ordering within each module at a critical temperature, followed 
by (ii) global ordering at a lower temperature. This indicates polarization of society into groups 
having contrary opinions can persist indefinitely even when mutual interactions between agents favor 
consensus. 
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Critical phenomena associated with order-disorder 
transitions is of central importance in statistical 
physics 1]. These results also have significant implica- 
tions for understanding social phenomena where coordi- 
nation dynamics is observed, e.g., in consensus forma- 
tion [2|, |3| and adoption of innovations Such pro- 
cesses can be analyzed using the spin models of statisti- 
cal physics which are generic systems for studying coop- 
erative phenomena. The spin orientations correspond to 
several equivalent but mutually exclusive choices made 
by an agent on the basis of information about the choice 
of the majority in its local neighborhood. The simplest 
case is when an agent decides between two competing 
choices where the dynamics can be modeled by an Ising 
system denned on a network reflecting the social contact 
pattern. 

One of the most remarkable features of complex net- 
works seen in many different contexts is their modular 
organization ||. Modular networks consist of subnet- 
works whose nodes have a significantly higher connection 
density compared to the overall density of the network 
S @| • A recent analysis of a network of mobile phone 
users, reconstructed according to their call frequency 
and duration, have shown the existence of such modular 
structure in the organization of social contacts . This 
has also been observed in other social networks, e.g., of 
scientific collaborators [icl |. electronic mail communica- 
tion [HG3, the PRETTY-GOOD-PRIVACY(PGP) en- 
cryption "web-of-trust" [l3| and even that of non-human 
animals |14|. Study of dynamical processes on such net- 
works |lfj | can help in understanding how individual be- 
havior at the microscopic level relates to social phenom- 
ena at the macroscopic level. 

In this Rapid Communication, we investigate a modu- 
lar network of Ising spins with ferromagnetic interactions, 
and report the existence of a phase with ordering among 
spins in each module, in the absence of global ordering. 
This modular ordered phase (Fig. [TJ a) is reached from 
the disordered state of the system through a continu- 
ous transition by lowering the temperature to the criti- 
cal value T™. Further reducing the temperature to Tj? 
results in another continuous transition to a state where 



all the modules are aligned in parallel, i.e., the globally 
ordered phase (Fig. [TJ b) . Our results help to understand 
how contrary opinions can co-exist in society even when 
mutual interactions favor consensus. In the social con- 
text, the concept of physical temperature corresponds to 
a measure of noise, e.g., due to imperfect information 
or uncertainty on the part of the agent, or the existence 
of idiosyncratic beliefs. Below the critical temperature 
Xjf, the relaxation time to global order diverges as the 
modularity increases. Thus, even when global order is 
achievable (i.e., T < T c 9 ), the time required to achieve 
consensus increases rapidly as the social organization be- 
comes more modular. However, this time can be altered 
significantly by considering processes that are observed 
in society. We implement this by (i) introducing an ex- 
ternal field modeling positive feedback effects which re- 
inforce the choice adopted by the majority or (ii) 
varying the strength of couplings between communities 
relative to those of intra-modular links, which is consis- 
tent with the well-known "weak tie" hypothesis for social 
networks [l7| . 

We consider a modular random network consisting 
of N nodes arranged into n m modules, each having n 
(= N/n m ) nodes [8]. The connection probability between 
a pair of nodes belonging to the same module is pi, while 
that for nodes belonging to different modules is p . The 
modularity of the network can be changed continuously 
by varying the ratio of inter-modular to intra-modular 
connectivity, r = ^ £ [0,1], keeping the average degree 
(k) constant. For r = 0, the network gets fragmented into 
a set of isolated clusters, while at r = 1, it becomes a ho- 
mogeneous or Erdos-Renyi random network. On placing 
Ising spins on the nodes of the network, the Hamiltonian 
for the system is given by 



H = 



JijO~i&j 

i,3 



(1) 



where a% = ±1 is the Ising spin on the ith node and the 
coupling strength Jy is J(> 0) if i,j are connected, and 
zero otherwise. The positive value of J (ferromagnetic 
coupling) implies that each link tries to align the two 
spins connected by it. 
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FIG. 1: (Color online) Ordering in modular networks: (a) 
Local ordering within modules without global order, and (b) 
the globally ordered state. 

Starting from a disordered state, we have performed 
Monte Carlo simulation for updating the spin states us- 
ing the Wolff algorithm 18]. For any configuration, wc 
can measure the average magnetic moment per module, 
fi — (| J2i£ S a i\)i t ne averaging being over all modules 
,s = 1, . . . , n m , and, the total (or global) magnetization 
of the network, M g = X^ '*- At equilibrium, the re- 
sulting phase diagram (Fig. [2]) clearly shows the exis- 
tence of three phases, corresponding to (a) the disor- 
dered state (n = 0, M g = 0), (b) a globally ordered state 
(fi ^= 0, M g ^ 0) and (c) a state with only modular or- 
der (p, 7^ 0, M g = 0). At low r, as the temperature is 
decreased the system undergoes two transitions: first, at 
T = T™, from the disordered state to the state with 
modular order, and next, at T = Tjf, to the globally or- 
dered state, where all spins in the network are aligned 
in parallel. As the modularity of the network decreases 
with increasing r, the two critical temperatures T™ and 
Xj? approach each other and finally converge. 

The phase diagram can be reproduced analytically by 
considering the free energy for the system. The magnetic 
moment of a single module is fj, = n(2f + — 1), where /+ 
is the fraction of "up" (a = +1) spins in that module. 
We assume that, at equilibrium, the magnetic moment 
of each of the n m modules has the same magnitude, with 
a fraction fT being positive (+fJL) and the remainder be- 
ing negative (— fj,). This is valid in the regime of strong 
modularity, i.e., r < 1. Thus, the total magnetic mo- 
ment for the system is M g = n m fi(2f" 1 — 1). The con- 
tribution to the internal energy of the system from each 
module is Ui — — JX,(2/ + — l) 2 , where Lj = "^" 2 ~ 1 ^ p i is 
the number of links within a module. This is based on 
the mean-field assumption that the neighborhood of all 
spins are identical. To obtain the internal energy con- 
tribution for interactions between modules, U , we note 
that each of the modules can be treated as "spins" of 
moment /i with L D = ""'("j" -1 ) p o links between them. 
Analogous to the preceding expression for Ui, we get 
U = -JL n 2 {2f™ - l) 2 . Then, the free energy for the 
system is 

F(f+J+) = n m (Ui - TSJ + U a - TS , (2) 
where the entropy terms, Si = — log (/_)_) + (1 — 




FIG. 2: (Color online) Critical behavior of Ising model on 
modular network, (a) The r-ksT/J phase diagram indicat- 
ing the existence of three distinct phases: disordered (M g = 
fj, = 0), modular order (M g = 0, (i / 0) and global order 
(M g ^ 0, /i / 0). The solid lines are phase boundaries ob- 
tained from Monte Carlo simulation, while the broken lines 
indicate the critical temperatures T™ and T| obtained from 
analytical expressions given in the text, (b) The distributions 
for global (m 9 = M g /N) and modular (m = fi/n) magnetiza- 
tions per site at different temperatures showing the two phase 
transitions for r = 0.002. At low temperature (T = 7J/fcg) 
both curves have bimodal nature, indicating global order. As 
T increases to lOJ/fcs (> T-!), the global order disappears 
resulting in unimodal distribution of m g . For temperatures 
higher than TJ™ (T = 15J/fcs), both curves are unimodal as 
the system is disordered. All results shown are for N = 512, 
n m = 16 and (k) = 14. 

/+)log(l - /+)], and S = -t^bL/T l°g(/+ ) + (1 - 
/™) log(l — /+*)], correspond respectively to the different 
ways in which nf + up spins can be distributed among 
n spins within each module, and n m /™ modules with 
moment +n can be distributed among n m modules. 

To analyze the critical behavior of the system we mini- 
mize the free energy F (Eq.[5J) with respect to /+, giving: 

1 , f+ T™ J Po n(n m -l)(2/;'~l) 2 

Insr = ■ 

4/+ - 2 1 — /+ T k B 

(3) 

where, T™ = 2 JLi/ (nks)- This indicates a continuous 
phase transition at the modular critical temperature TJ", 
below which spins within a module are ordered but f" 1 = 
1/2. This state corresponds to the phase for which // ^ 
0,M g = 0. As the temperature decreases below T™, 
the different modules get aligned with each other at a 
temperature Tf. This is obtained by minimizing F with 
respect to /™: 

_ W — i± — = 1°. (4) 

2(2/™ - 1) 6 1 - f™ T ' w 

which shows a continuous phase transition at the global 
critical temperature Tf> = f^(n m - l)n 2 (2/ + - l) 2 . The 
expression for Tf! does not liave a closed analytic form 
and it is obtained by numerical minimization. As r — > 1, 
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the network loses its modular structure and becomes a 
homogeneous random network, so that Eqs. OH] give the 
critical temperature T c m = T c 9 = J{k)/k B - 

The above analysis for finite-size networks are valid 
even in the thermodynamic limit when it is approached 
by increasing the number of modules, n m . Note that, if 
we instead increased the number of nodes in a module, 
n (keeping n m fixed), the modularity of the network is 
lost and hence, there is only a single continuous order- 
disorder transition at T c = J{k)/kB- The two types of 
ordering seen for the Ising model on a modular network 
imply that not only can consensus formation in a society 
be an extremely slow process , but under certain con- 
ditions it may never be achieved and communities with 
contrary opinions can persist indefinitely. 

We now look at how the time required to reach the 
equilibrium state corresponding to global order varies 
with the modularity of the network. Starting from an 
initially disordered state (M g = 0, A t = 0), spins are up- 
dated using Glauber dynamics (201 . Fig. [3J shows that 
the relaxation to global order takes an extremely long 
time as the system become more modular. To under- 
stand this, we first consider a single isolated module of n 
nodes with Li intra- modular links. Under the mean- field 
approximation, the time-evolution of the magnetization 
per site m = fi/n is described by 



dm(t) 
dt 



+ m(t) — tanh 



(5) 



where T c = 2JLi/(nks)- At equilibrium, this reduces 
to the usual mean-field equation, m = tanh(mT c /T), in- 
dicating that T c is the critical temperature. Defining 
relaxation time r to be the time in which m becomes 0.5, 
we obtain from Eq. \5\ 



0.5 



dx 



o+ 



tanh(a:T c /T) 



(6) 



Although it is small at low temperatures, r diverges as 
T — > T c due to critical slowing down. 

In a network consisting of several modules, the relax- 
ation time to modular order, r mm , is identical to r de- 
rived above with T c = T c m . To obtain the relaxation 
time to global order, we first note that the free energy of 
a single module F m (f + ) — Ui — TSi, has two minima at 
/+ and 1 — /+ (say). These correspond to the magnetic 
moments ±/i, which are separated by a free energy max- 
imum at f + = 1/2 corresponding to fi = 0. To switch 
from +/i to — n (or vice versa), a module has to over- 
come an energy barrier A = F m {l/2) — F m (f+). Thus, 
the time required to attain a global magnetization M g is 
slowed down by the factor exp[A/(fcsT)]. Defining the 
global relaxation time, r gm , to be the time in which the 
global magnetization per site, M g /(n m fi) becomes 0.5, 
we obtain from Eq. [6l 



exp {if) 



0.5/|2/°-l| 



dx 



o+ 



tanh(xT c 9 /r) 



(7) 
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FIG. 3: Relaxation time to the globally ordered state, r 9m , 
as a function of r at different temperatures. With increasing 
modularity, r gm diverges. (Inset) The variation with 
1/T at r = 0.008. For all cases, N = 512, n m = 16 and 
(ft) = 14. 



Here we have assumed that T gm 3> t„ 



which is valid 



for T ^ T 9 and low r. In this region, r gm diverges with 
increasing modularity of the network, the trend agree- 
ing with the numerical results of Fig. [3] Note that, at 
low temperatures (T <C Tg), the integral in Eq. [7] can 
be neglected and r gm ps exp(A//csT) cx exp(Jpi/fcs? 1 ), 
assuming « 1. Thus, the relaxation to global or- 
der becomes exponentially slower as the temperature de- 
creases (Fig. [3j inset). The above results indicate that 
even when T < T 9 , a strongly modular network may re- 
quire an extremely long time to reach the globally ordered 
state, which explains why previous studies may have er- 
roneously observed a separate phase without global order 
in this region [2l| . 

It may appear from the preceding analysis that achiev- 
ing global consensus is extremely difficult in a real social 
network having modular organization. However, we now 
discuss possible mechanisms by which the time to attain 
global order can be changed significantly. First, we look 
at the role of an external magnetic field which is propor- 
tional to the global magnetization M g . This corresponds 
to positive feedback effects in social systems, where, al- 
though two competing alternatives are initially equiva- 
lent, as more and more agents switch to one particular 
option, it becomes the preferred choice [l6j]. Introduc- 
ing such a field, the Hamiltonian for the system becomes 
H = Sij JijO~iCfj — ^(53i cr i) 2 - The external field h has 
no effect in the absence of global order, but when M g is 
non-zero, the field drives the system to the equilibrium 
state corresponding to global order much faster as seen in 
Fig. 2](a). The free energy in this case can be obtained by 
replacing JL m with JL m + hn 2 m in Eq.[5] Thus, the field 
effectively increases the inter-modular interactions (mak- 
ing the network less modular) which drives the system 
away from the critical point by increasing Xj?, thereby 
reducing T gm . 

We next consider the possibility that the interaction 
strengths for inter-modular connections ( J ) may be dif- 
ferent from those of intra- modular links (Ji). From 
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FIG. 4: The relaxation time to the globally ordered state, r 9m , 
as a function of (a) the external field h, and, (b) the ratio of 
inter to intra-modular coupling strength J /Ji, at different 
temperatures. The minima of T gm occur when Jo/Ji — (kin), 
the average number of links a node has with other nodes in 
its module. For all cases, JV = 512, n m = 16, (k) — 14 and 
r = 0.004. 



Eq. [2l it is clear that increasing the ratio J /Ji is equiv- 
alent to increasing r. Therefore, as the intra-modular 
links increase in strength relative to the inter-modular 
links (making the network even more modular), the time 
to achieve global order (r gm ) increases. On the other 
hand, when the inter-modular links are stronger, T gm 
decreases up to a point and then, with increasing J Q , 
starts increasing as the nodes in different communities 
get strongly coupled thereby destroying the identity of 
individual modules (Fig. |4j b). This non-monotonic be- 
havior of relaxation time as a function of the ratio of 
strengths for short- and long-range interactions is in con- 
trast with that seen in the case of Watts-Strogatz (WS) 
small- world network model [12] • This reinforces previous 
results that the dynamics of WS networks are strikingly 



different from that of modular networks, although their 
structural properties are similar H . It also suggests that 
a system may maintain diversity by using weak links be- 
tween their constituent communities 

In this Rapid Communication we have shown that the 
order-disorder transition in an Ising model defined on 
a modular network shows the existence of three distinct 
phases. While the disordered and globally ordered phases 
are similar to those expected for a homogeneous network, 
the existence of a phase with local ordering within mod- 
ules but no global order is a novel effect of the modular 
structure. This has significant implications for ordering 
dynamics in real systems, e.g., consensus formation on 
social networks. It indicates that, under certain condi- 
tions, homogeneous groups with contrary opinions can 
coexist forever even when mutual interactions between 
agents favor consensus, so that society becomes polar- 
ized. Our results suggest that such tendencies can be 
countered by increasing inter-community communication 
and improving the overall penetration of mass media in 
the society. Although the present Rapid Communication 
looks at the case of two competing choices, it is possible 
that the effects seen here extend to the case of agents 
choosing between multiple (> 2) alternatives, e.g., in the 
context of q— state Potts spin dynamics on modular net- 
works. Further, the network can be made more realistic 
from a social perspective by considering the presence of 
hubs and a broad link-strength distribution Q. In sum- 
mary, our results imply that due to mesoscopic structural 
inhomogeneity, equilibrium as well as dynamical proper- 
ties of local regions in a complex network may depart 
significantly from those for the entire network. 
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